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Executive  Summary 


This  paper  describes  the  procedures  followed  to  develop  a  detailed  watershed 
database  for  Fort  Benning,  Georgia,  and  the  surrounding  area.  In  addition  the 
relationship  between  watershed  morphology  and  stream  order  was  examined. 
Watershed  order  and  a  number  of  variables  describing  surface  topography  and  the 
stream  network  were  computed  and  statistical  analysis  procedures  were  used  to 
develop  a  predictive  relationship.  Watershed  boimdaries  were  computed  from  a 
digital  elevation  model  and  assigned  an  order  using  the  Strahler  stream-ordering 
technique.  This  procedure  is  rather  tedious  and  requires  that  all  the  drainage 
network  upstream  of  the  area  of  interest  be  digitized.  A  number  of  physical 
parameters  defining  these  watersheds  were  computed  using  a  geographic 
information  system  (GIS)  and  the  relationship  of  these  parameters  to  stream  order 
was  examined.  The  purpose  of  this  analysis  was  to  determine  if  stream  order  could 
be  predicted  reliably  using  a  number  of  computed  physical  parameters.  Regression 
analysis  showed  that  stream  order  had  a  highly  significant  relationship  (r-square  = 
0.77)  with  two  easily  computed  variables:  total  rehef,  and  average  slope.  A 
procedure  for  estimating  stream  order,  within  an  acceptable  degree  of  error,  could 
be  beneficial  for  many  applications,  such  as  assisting  in  the  parameterization  of 
hydrologic  models. 


1  Introduction 


Background 

The  Department  of  Defense  (DoD)  has  established  ecosystem  management  as 
its  approach  to  managing  military  lands.  The  Strategic  Environmental  Research 
and  Development  Program  (SERDP),  Ecosystem  Management  Program  (SEMP), 
was  established  in  December  1997  to  help  address  critical  deficiencies  in 
knowledge  that  prohibit  the  DoD  from  fully  achieving  this  goal  (SERDP  2000). 
One  component  of  the  SEMP  is  the  Ecosystem  Characterization  and  Monitoring 
Initiative  (ECMI).  The  objectives  of  the  ECMI,  described  in  SERDP  (2000), 
include  developing  and  implementing  methods  of  monitoring  and  characterizing 
ecosystems  that  can  help  land  managers  assess  relationships  between  land  use  and 
management,  and  ecosystem  structure,  function,  and  pattern.  Fort  Benning, 
Georgia,  was  established  as  the  first  test  site. 

Watersheds,  which  provide  linked  gradients  of  terrestrial,  riparian,  and 
aquatic  systems,  were  established  as  the  critical,  ecosystem  delineation/m^ping 
xmit.  A  watershed  (also  known  as  a  hydrologic  unit,  catchment,  or  drainage  basin) 
is  defined  as  that  area  of  land  draining  into  a  particular  stream  or  other  surface 
water  body.  For  any  location  in  a  stream,  there  is  an  associated  area  that 
contributes  water  to  its  flow.  The  watershed  divide  is  that  line  which  divides  the 
area  contributing  water  to  the  stream  and  that  which  contributes  water  to 
neighboring  streams  or  water  bodies.  Therefore,  each  watershed  is  defined  by  its 
outlet  or  pour  point  (the  point  in  the  stream  which  receives  all  water  in  the 
watershed)  and  the  associated  watershed  divide  derived  from  that  point  and  the 
local  topography. 

The  primary  goal  of  this  project  was  to  develop  a  detailed  digital  elevation 
model  (DEM)  and  watershed  boundaries  to  support  the  efforts  of  the  SEMP 
ECMI.  A  secondary  goal  was  to  examine  the  relationship  between  stream  order 
and  watershed  geomorphology.  This  paper  documents  the  process  by  which 
watersheds  were  produced  for  the  Fort  Benning  study  area,  as  well  as  a  statistical 
analysis  of  the  relationship  between  stream  order  and  watershed  characteristics. 

Although  a  number  of  early  studies,  such  as  Morisawa  (1959),  have  analyzed 
the  relationship  between  quantitative  geomorphology  and  stream  measurements, 
most  of  these  studies  were  limited  by  the  complexity  of  calculating  detailed 
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watershed  variables.  The  development  of  geographic  information  systems  (GIS) 
has  made  it  possible  to  quickly  define  watershed  boundaries  and  to  compute 
geomorphic  variables  describing  them. 


Study  Site  Description 

Fort  Beiming  is  located  in  west-central  Georgia,  south  of  the  city  of 
Columbus,  Georgia.  It  occupies  approximately  73,533  hectares  (ha)  in 
Chattahoochie,  Muscogee,  and  Marion  Counties  in  Georgia,  and  Russell  County 
in  Alabama.  The  base  lies  within  the  humid  temperate  domain,  subtropical 
division,  coniferous-broadleaved  semi-evergreen  forest  province,  as  defined  by 
Bailey  (1995).  Fort  Benning  falls  within  the  southeastern  plains  ecoregion  as 
defined  by  Omemik  (1987).  The  base  is  located  within  hydrologic  unit 
03 130003,  as  defined  by  the  U.S.  Geological  Survey  (USGS)  (Figure  1). 

The  Fort  Benning  SEMP  demonstration  site  is  described  in  much  greater 
detail  in  “Design  Document  for  Long-Term  Monitoring  Program,  Fort  Benning, 
Georgia”  (SERDP  2000). 
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Figure  1 .  Fort  Benning  and  associated  hydrologic  unit 
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2  Database  Development 


The  development  of  watershed  boundaries  requires  two  types  of  data:  detailed 
topographic  information,  and  an  accurate  and  topologically  correct  stream  net¬ 
work.  For  this  study,  both  types  of  data  were  derived  from  USGS  1 :24,000  maps. 

Arc/Info  GIS  software  (ESRI 2000)  was  used  to  develop  the  database  and  for 
subsequent  data  processing. 


Elevation  Contours  and  Benchmarks 

Digital  versions  of  vector  contour  lines  (Figure  2)  were  purchased  from 
Landinfo  International  for  each  of  the  1: 24,000-scale  USGS  topographic 
quadrangle  maps  covering  the  hydrologic  unit  in  which  Fort  Benning  is  situated. 

The  contour  lines  were  appended  into  one  file  and  the  following  steps  were 
taken  to  ensure  that  the  data  were  accurate: 


a.  The  elevations  of  the  contours  lines  were  checked  closely  to  verify  that 
they  represented  the  correct  range  of  values  for  each  source  map  sheet. 

b.  The  digital  files  were  plotted  on  clear  Mylar  overlays  and  compared  with 
the  paper  source  maps. 

c.  The  files  were  inspected  for  contour  lines  carrying  incorrect  elevation 
attributes. 

d.  Contours  that  were  located  along  edges  of  the  source  maps  were 
edgematched  and  node  errors  were  corrected. 


To  supplement  the  contour  information,  elevation  benchmarks  were  digitized 
from  each  of  the  USGS  quadrangle  maps  and  the  proper  elevations  were  assigned 
as  attributes  to  the  resulting  point  coverage. 


Stream  Network 

For  this  project,  streams  were  identified  as  blue  lines  on  the  USGS  1:24,000 
scale  maps,  as  required  by  NI  170-304  SubPart  C,  304.20(a)  (NRCS,  1995).  The 
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Figure  2.  Vector  topographic  contours  over  portion  of  Fort  Benning,  Georgia 


smallest  streams  were  those  with  no  tributaries,  while  the  largest  streams 
discharged  directly  to  the  Chattahoochie  River. 

A  tablet  digitizer  was  used  to  obtain  the  streams  from  each  USGS  map.  The 
streams  were  then  edgematched  to  ensure  that  they  were  properly  connected. 
Check  plots  were  produced  to  compare  the  digital  versions  of  the  stream  segments 
to  the  source  1:24,000  maps. 

The  following  steps  were  then  followed  to  prepare  the  streams  for  the 
ordering  process  and  Ae  generation  of  an  elevation  model. 
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a.  All  stream  segments  were  checked  to  ensure  that  they  were  connected  and 
oriented  (pointing)  in  a  downstream  direction.  This  is  required  for  the 
stream  ordering  procedure. 

h.  The  files  were  checked  for  accurate  topology  and  to  ensure  that  there 
were  no  overlapping  line  segments. 

c.  The  streams  were  also  inspected  closely  to  assure  that  they  crossed 
contours  at  the  appropriate  locations. 

d.  The  stream  segments  were  densified,  i.e.,  the  number  of  vertices  in  each 
stream  segment  was  increased. 

e.  In  areas  where  lakes  or  small  ponds  interrupted  stream  segments  on  the 
source  maps,  line  segments  had  to  be  created  to  join  both  the  inlet  and 
outlet  streams  into  a  continuous  line.  This  is  a  required  for  the  stream 
ordering  procedure. 
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Elevation  Model  Development 

The  TOPOGRID  software  in  Arc/lnfo  (ESRI,  2000)  was  used  to  produce  the 
final  surface  elevation  model.  TOPOGRID  is  an  interpolation  method, 
specifically  designed  to  generate  hydrologically  correct  OEMs  from  elevation  data 
and  stream  networks.  It  is  based  on  the  ANUDEM  program  developed  by 
Michael  Hutchinson  (1988, 1989). 


The  data  used  to  develop  the  elevation  model  are  listed  in  Table  1; 


Table  1 

Data  Used  to  Create  the  Digital  Elevation  Model 

Coverage  Name 

Topology 

Description 

Contours 

Line 

Elevation  contours  from  USGS  1 :24,000  maps 

Benchmarks 

Point 

Elevation  benchmarks  from  USGS  1 :24,000  maps 

Streams 

Line 

Topologically  correct  stream  network 

Boundary 

Polygon 

Polygon  representing  boundary  of  interpolation 

The  following  parameter  file  was  used  to  generate  the  elevation  model: 

TOPOGRID  output  surface  10 
ENFORCE  ON 
DATATYPE  CONTOUR 
MARGIN  0.0 
ITERATIONS  30 
TOLERANCES  2.5  1.0  0.0 
XYZLIMITS##  #### 

CONTOUR  contours  ELEV-M 
POINT  benchmarks  ELEV-M 
STREAM  streams 
BOUNDARY  boundary 
OUTPUTS  sinks  drainage 
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A  shaded  relief  depiction  of  the  final  elevation  model  is  presented  in  Figure  3  and 
a  portion  of  the  elevation  model  is  shown  in  a  ‘‘hill-shaded”  format  in  Figure  4. 


Figure  3.  Shaded  relief  depiction  of  elevation  model 
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Figure  4.  Stream  network  (blue)  over  a  hillshade  depiction  of  the  digital  elevation  model 


Stream  Ordering  Procedure 

Stream  ordering  is  a  process  of  identifying  and  grouping  stream  segments  and 
their  corresponding  watersheds  in  terms  of  size  and  complexity.  Theoretically, 
watersheds  of  similar  order  display  similar  hydraulic  properties  and  ecological 
function.  There  are  four  commonly  described  approaches  to  stream  ordering.  In 
this  study,  the  approach  used  was  that  originally  described  by  Horton  (1945)  and 
revised  by  Strahler  (1952).  In  this  ordering  scheme,  the  smallest  stream  segments 
near  the  drainage  divide  are  assigned  the  lowest  order  (i.e.,  first-order  stream)  and 
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the  stream  segment  at  the  watershed  outlet  is  assigned  the  highest  order.  Each 
sub-basin  identified  is  assigned  the  same  order  as  the  largest  stream  segment 
within  it. 


The  ordering  system  can  be  described 
by  the  following  series  of  steps  and  is 
depicted  graphically  in  Figure  5  (Chow, 
Maidment,  and  Mays  1988): 

a.  The  smallest  recognizable 
channels  are  designated  order  1; 
these  channels  normally  flow 
only  during  periods  of  wet 
weather. 

b.  Where  two  channels  of  order  1 
join,  a  channel  of  order  2  results 
downstream;  in  general,  where 
two  channels  of  order  /  join,  a 
channel  of  order  /  +  1  results. 

c.  Where  a  channel  of  lower  order 
joins  a  channel  of  higher  order, 
the  channel  downstream  retains 
the  higher  of  the  two  orders. 

d  The  order  of  the  drainage  basin  is  designated  as  the  order  of  the  stream 
draining  its  outlet,  the  highest  stream  order  in  the  basin,  I. 

Every  stream  on  the  map  could  be  ordered  as  long  as  its  furthest  upstream 
extent  was  known.  This  ordering  was  completed  for  all  streams  over  30  m  (100  ft) 
in  length  on  the  USGS  quadrangle  maps  on  the  Georgia  side  of  Fort  Benning,  but 
was  not  done  for  the  Chattahoochie  River  itself  (as  the  area  drained  by  this  major 
river  extends  beyond  the  hydrologic  unit  study  area). 

This  process  was  conducted  for  the  Fort  Benning  drainage  system  using  an 
automated  Arc  Macro  Language  (AML)  script  operating  within  the  Arc/Info  GIS 
software  package.  The  program  used  to  order  the  streams  is  listed  in  Appendix  A. 


Watershed  Boundary  Delineation 

Several  procedures  within  the  GIS  software  were  used  to  produce  polygons 
defining  the  drainage  area  associated  with  each  stream  segment  or  group  of  seg¬ 
ments.  These  polygons  represent  the  various  watersheds  having  different  orders. 
Several  steps  were  required  by  produce  these  watershed  boundaries  from  the 
elevation  model. 


Figure  5.  Diagram  of  the  Strahler  stream 
ordering  procedure 
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First,  the  Arc/Info  GRID  function  FLOWDIRECTION  was  used  to  generate  a 
raster  file  depicting  flow  direction  from  each  cell  in  the  elevation  model  to  its 
steepest  downslope  neighbor.  The  method  of  deriving  flow  direction  employed  by 
the  Arc/Info  software  is  described  in  Jenson  and  Dominique  (1988). 

The  watersheds,  as  described  previously,  are  partly  defined  by  their  outlets. 
Points  where  two  streams  intersect  define  watershed  outlets.  The  nodes  where 
streams  intersected  were  selected  from  the  stream  coverage  and  copied  into  a 
separate  Arc/Info  coverage.  This  coverage  of  outlets  was  then  converted  to  a  raster 
format. 

Finally,  the  watershed  boundaries  were  defined  using  the  Arc/Info  GRID 
WATERSHED  function,  with  both  the  output  of  the  FLOWDIRECTION  process 
and  the  raster  dataset  of  outlets  used  as  input.  The  resulting  watershed  grid  was 
then  converted  from  a  raster  file  into  a  polygon  coverage.  Manual  editing  was 
required  to  clean  up  the  final  watershed  polygons  to  remove  artifacts  such  as 
polygons  representing  very  small  areas,  i.e.,  one  cell  in  the  source  raster  dataset. 

The  next  step  was  to  assign  orders  to  the  watershed  polygons.  The  ordering  of 
watersheds  follows  the  ordering  of  the  streams  that  drain  them.  If  the  outlet  of  an 
order  i  stream  is  that  point  where  it  joins  another  stream  of  order  i  or  higher,  an 
order  i  watershed  corresponding  to  that  stream  can  be  derived  from  the  surround¬ 
ing  topography.  Although  each  reach  of  a  stream  will  have  a  unique  order,  each 
area  of  land  may  belong  to  more  than  one  order  of  watershed.  In  general,  any 
order  i  watershed  will  contain  at  least  two  watersheds  of  order  /  - 1,  and  each 
order  /  watershed  will  be  contained  in  some  watershed  of  higher  order.  The 
IDENTITY  command  assigned  the  stream  coverage  order  to  the  watershed 
polygons.  This  process  was  conducted  for  each  order  of  stream,  resulting  in  six 
different  polygon  coverages. 

A  total  of  3348  watershed  boundaries  were  delineated.  The  numbers  of 
watersheds  (by  order)  delineated  in  the  study  area  are  listed  in  Table  2.  Maps  of 
each  watershed  order  are  presented  in  Figures  6-11. 


Table  2 

Number  of  Watersheds  Delineated  in  the  Fort  Benning  Study  Area 
(by  order) 

Watershed  Order 

Number  of  Watersheds 

1 

2570 

2 

596 

3 

142 

4 

31 

5 

8 

6 

1 
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Figure  6.  First  order  watersheds 
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Figure  7.  Second  order  watersheds 
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Figure  8.  Third  order  watersheds 


Figure  9.  Fourth  order  watersheds 
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Figure  10.  Fifth  order  watersheds 
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4  Statistical  Analysis  of 
Watershed  Order 


The  goal  of  this  analysis  was  to  relate  a  number  of  physical  variables  to 
watershed  order.  These  variables  (Table  3)  are  fairly  straightforward  to  measure 
within  a  GIS.  The  parameters  that  were  measured  are  described  in  Table  3; 


Table  3 

Variables  Used  in  the  Statistical  Analysis 

Variable  (Variable  name) 

Type  of  Variable 

Unit  of  Measure 

Description 

Order  (ORDER) 

Interval 

_ 

Unitiess 

Dependent  Variable 

Area  (BASIN_AREA) 

Ratio 

Square  meters 

Total  area  of  drainage  area 
above  an  outlet 

Number  of  Stream  Segments 
(NUM_STREAM) 

Interval 

Integer  value 
(unitless) 

Count  of  total  number  of 
stream  segments  upstream 
from  a  watershed  outlet 

Total  Length  of  Stream 
Segments  (TOT_LENGTH) 

Ratio 

Meters 

Total  length  of  stream 
segments  upstream  from  a 
watershed  outlet 

Perimeter  (BASIN_PERIM) 

Ratio 

Meters 

Total  perimeter  of  drainage 
area  above  an  outlet 

Maximum  Relief 
(TOT^RELIEF) 

Ratio 

Meters 

Max  elevation  -  min 
elevation 

Average  Slope 
(AVG_SLOPE) 

Ratio 

Degrees  slope 

Average  slope  in  the  area 

Method 

Of  all  the  watersheds  listed  in  Table  1, 30  from  each  order  (orders  1  to  4) 
were  randomly  selected  from  the  GIS  database.  This  was  done  to  make  the 
sample  from  each  order  as  equal  as  possible  to  prevent  the  first  order  watersheds 
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from  biasing  the  analj^is.  The  random  selection  process  resulted  in  a  sample  of 
129  watersheds  for  the  regression  analysis.  Variables  listed  in  Table  3  were 
calculated  using  the  Arc/Info  GIS  software  (ESRI 2000).  The  complete  dataset  is 
attached  (Appendix  B). 

A  multiple  regression  approach,  with  stepwise  backwards  variable  selection, 
was  used  to  determine  which  independent  variables  were  significant  in 
determining  stream  order. 


Results  and  Analysis 


Model  1 

In  a  model  that  included  all  six  computed  variables,  four  were  foimd  to 
contribute  significantly  to  the  first  model  (basin  area,  basin  perimeter,  total  relief, 
and  average  slope).  The  model  resulted  in  an  r-square  of 0.7741.  It  was 
somewhat  surprising  that  neither  of  the  variables  directly  related  to  the  stream 
network  (number  of  streams  and  total  length  of  streams)  was  found  to 
significantly  add  to  the  model. 


Table  4 

Model  1  Results 

Variable 

Parameter 

Estimate 

Standard  Error 

t  value 

Pr>|t| 

Variance 

Inflation 

Basin  Area 

-6.87469  X  10  ® 

1.558417x10-® 

-4.41 

<  0.0001 

9.60798 

Basin  Perim 

0.00003411 

0.00000727 

4.69 

<  0.0001 

19.88638 

Total  Relief 

0.01528 

0.00356 

4.30 

<  0.0001 

5.89787 

Avg  Slope 

-0.10353 

0.04209 

-2.46 

0.0153 

1.21827 

Close  examination  of  this  complete  model  shows  a  significant  problem  with 
multicolinearity.  Variance  inflation  values  are  very  high,  indicating  that  several  of 
the  variables  are  highly  correlated.  Eigenvalues  and  condition  index  values  are 
shown  in  Table  5. 


Table  5 

Model  1  Colinearity  Diagnostics 

Variable 

Eigenvalue 

Condition  Index 

Basin  Area 

1.18555 

1.75120 

Basin  Perim 

0.13026 

5.28317 

Total  Relief 

0.03451 

10.26369 

Avg  Slope 

0.01395 

16.14182 
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Model  2 


Given  the  amount  of  colinearity  present  in  Model  1,  a  number  of  alternate 
models  were  tested  by  removing  variables.  This  process  of  elimination  showed 
that  a  model  that  included  only  two  of  the  variables  (total  relief,  and  average 
slope)  (Table  6)  accounted  for  almost  all  the  total  variation  of  Model  1  (r-square 
of  0.7339),  while  eliminating  colinearity  problems  (Table  7). 


Table  6 

Model  2  Results 

Parameter 

Variance 

Variable 

Estimate 

Standard  Error 

t  value 

Pf>i(i 

inflation 

Total  Relief 

0.02969 

0.00160 

18.59 

<  0.0001 

1.02589 

Avg  Slope 

-0.17792 

0.04159 

-4.28 

<  0.0001 

1.02589 

Table  7 

Model  2  Colineaiity  Diagnostics 

Variable 

Eigenvalue 

Condition  Index 

Total  Relief 

0.11886 

4.89377 

Avg  Slope 

0.03461 

9.06855 

Examination  of  observation  diagnostics  (DFFITS  and  DFBETAS)  revealed 
that  a  number  of  observations  were  anomalous.  The  critical  value  of  DFFITS 
(2*sqrt(p/n)  =  0.25)  and  DFBETAS  (2/(sqrt(«))=0.176)  revealed  a  number  of 
irregular  observations.  Also,  the  RSTUDENT  diagnostic  was  used  to  evaluate 
outliers.  Results  of  these  analyses  are  presented  in  Table  8. 

Obviously,  observations  should  not  be  deleted  randomly  to  improve  the 
performance  of  a  model.  However,  each  of  the  questionable  basins  was  displayed 
in  the  GIS,  which  showed  that  observations  18  and  103  were  located  in  the 
easternmost  portion  of  the  study  area.  The  streams  in  this  area  were  digitized 
from  reference  maps  that  differed  somewhat  from  the  rest  of  the  study  area.  These 
maps  were  produced  during  a  different  period  from  the  rest  of  the  maps  and  the 
streams  on  these  source  maps  were  noticeably  less  detailed  than  those  on  the 
source  maps  that  had  been  used  throughout  the  rest  of  the  study  area.  Therefore, 
these  two  observations  were  deleted  from  the  model. 
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Table  8 

Influential  Observations  and  Outliers  (numbers  given  represent 
observation  numbers  in  data  listing  -  Appendix  B) 

DFFITS 

DFBETAS 

RSTUDENT 

crit  vaf  =  0.25 

crit  val  =  0.176 

crit  val  =  2.7 

4 

18 

18 

25 

25 

43 

43 

88 

88 

101 

101' 

103 

103 

103’ 

115 

116 

127 

^  Observations  very  close  to  but  not  greater  than  critical  value. 

Model  3 

Tables  9  and  10  (the  final  model)  show  the  results  of  eliminating  the  two 
questionable  observations.  The  final  model  produced  an  r-square  of 0.7684, 
which  is  nearly  as  good  as  Model  1,  which  had  four  variables.  There  is  no 
problem  with  multicolinearity  in  this  model  and  observation  and  outlier 
diagnostics  look  very  good. 

The  SAS  program  is  presented  in  Appendix  C. 


Table  9 

Model  3  Results 

Parameter 

Variance 

Variable 

Estimate 

Standard  Error 

t  value 

Pr>\tl 

Inflation 

Total  Relief 

0.02998 

0.00149 

20.17 

<  0.0001 

1.02333 

Avg  Slope 

-0.16417 

0.03909 

-4.20 

<  0.0001 

1,02333 

Table  10 

Model  3  Colinearity  Diagnostics  | 

Variable 

Eigenvalue 

Condition  Index 

Total  Relief 

0.12050 

4.85925 

Avg  Slope 

0.03417 

9.12470 
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A  Shapiro-Wilk  analysis  was  conducted  to  examine  the  residuals  for 
normality.  Results  are  shown  in  Table  1 1. 


Table  11 

Model  3  Residual  Normality  Test 

Test 

Statistic 

p  Value 

Shapiro-Wilk 

W  0.99492 

Pr<W  0.9332 

Table  12  summarizes  the  results  of  each  of  the  models  presented  above. 


Table  12 

Results  of  Stepwise  Backwards  Multiple  Regression  Analysis 

Model 

(Case) 

Regression  Equation 

Std.  Error  of  Y 

1 

-6.87(BA)  +  0.00003411  (BP)  +  0.01528(TR) 

+  0.01395  (AS) 

0.7741 

0.22912 

2 

0.02969(TR)  -  0.17792(AS) 

0.7339 

0.24058 

3 

0.02998(TR)  -  0.1 641 7(AS) 

0.7684 

0.22764 

Note:  BA 

=  Basin  Area;  BP  =  Basin  Perimeter;  TR  =  Total  Relief;  AS  =  Average  Slope.  | 
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5  Conclusions 


The  primary  objectives  of  this  project  were  to  develop  a  high-resolution  DEM 
and  to  produce  watershed  boimdaries  and  an  ordered  stream  network  from  this 
DEM.  As  the  SERDP  ECMI  monitoring  plan  is  watershed-based,  these  products 
were  critical  to  support  both  monitoring  design  efforts  and  fiiture  research  efforts. 
The  elevation  model,  stream  network,  and  watershed  boundaries  have  been  placed 
on  the  ECMI  Internet-based  data  repository  and  are  available  to  SEMP 
researchers. 

A  secondary  goal  of  this  project  was  to  statistically  analyze  of  the  watersheds 
to  gain  a  better  understanding  of  the  geomorphic  relationships  between  stream 
order  and  associated  drainage  areas.  Strong  relationships  were  developed  between 
total  relief  and  average  slope. 

Surprisingly,  a  model  that  included  only  total  relief  and  average  slope  gave 
the  best  results.  Before  this  analysis,  basin  area,  the  number  of  stream  segments, 
and  the  total  length  of  streams  were  thought  to  be  much  more  important  than  relief 
and  slope.  Interestingly,  this  model  requires  absolutely  no  information 
whatsoever  regarding  the  stream  network.  In  addition,  it  is  very  easy  to  compute 
slope  and  relief  using  a  GIS.  Therefore,  this  model  could  prove  useful  for  quickly 
determining  the  order  of  a  stream  in  a  study  area,  without  requiring  the 
compilation  of  detailed  stream  data  further  up  the  drainage  basin. 
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Appendix  A 

Stream  Ordering  Program 


/*  This  program  riins  from  within  Arc/Info,  specifically  within 
/*  the  Arcedit  module  of  Arc/ Info. 

/*  It  is  written  in  Arc  Macro  Language 
/* 

&args  cover 

/*  shreve  and  strahler  ordering  a  coverage 

&if  [show  program]  ne  * ARCEDIT*  &then  &do 
&type  this  starts  from  arcedit 
&pause  &seconds  3 
areturn 
&end 

/*  assume  the  item  JORDER  numbers  the  arcs  from  the  outlet, 

/*  and  that  there  are  no  pseudo-nodes 
ec  % cover % 
ef  arc 

Seif  [iteminfo  %cover%  -ARC  JORDER  -exists]  &then 

&return  *This  program  expects  %cover%  to  have  the  item  JORDER 
&if  [iteminfo  %cover%  -ARC  SHREVE  -exists]  &then 
additem  SHREVE  441 

&:if  [iteminfo  %cover%  -ARC  STRAHLER  -exists]  &then 
additem  STRAHLER  221 
sel  all 

calc  %cover%-ID  =  %cover%# 
calc  SHREVE  =  0 
calc  STRAHLER  =  0 

/*  statistics 
/*  maximum  JORDER 

/*  maximum  %cover%-ID 

/*  END 

/*  aset  ord  [show  statistic  1  1] 

/*  &set  maxid  [show  statistic  2  1] 

/*  &type  %maxid%  arcs  up  %ord%  levels  from  the  mouth. 

sel  dangle 
resel  j order  ne  1 
calc  SHREVE  =  1 

calc  STRAHLER  =1  /*  all  sources  are  coded 

sel  SHREVE  =  0 
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&label  bigloop  /*  look  at  every  imcoded  arc 
sel  SHREVE  =  0 
&s  num  [show  number  select] 

&ty  We  now  have  %num%  uncoded  arcs 
&if  %num%  eq  0  &then 
&goto  done 

&do  index  =  1  &:to  %num% 

ScS  id  =  [show  arc  [show  select  %index%]  item  %cover%-ID] 

&set  id%index%  -  %id%  /*  write  each  one  to  an  array 
&end 

ado  index  =  1  &to  %num%  /*  for  each  uncoded  arc 
&set  id  =  [value  id%index%] 
sel  %cover%-ID  =  %id% 

/*  &ty  ID  %index%  is  %id% 

&s  jo  =  [show  arc  [show  select  1]  item  jorder] 
select  connect 
resel  jorder  gt  %jo% 

&s  upnum  =  [show  number  select] 

&:if  %upnum%  eq  0  &then 

Sere  turn  Error  looking  upstream  from  %cover%-ID,  arc  %id% 
resel  shreve  gt  0 

Self  [show  number  select]  It  %upnum%  &then 
Segoto  continue 
&s  this shreve  =  0 
&:S  maxstrahler  =  0 

/*  we  can  handle  more  than  two  upstream  arcs,  though  results  may  not  be 
perfect 

&:S  ui  =  1  /*  build  my  own  loop 

Silabel  innerloop  /*****&do  ui  =  1  &to  %upnum% 

ScS  ushreve  =  [show  arc  [show  select  %ui%]  item  shreve] 

/*  &:if  %ushreve%  eq  0  athen; Siigoto  continue 

ScS  thisshreve  =  %thisshreve%  +  %ushreve% 

ScS  ustrahler  =  [show  arc  [show  select  %ui%]  item  strahler] 

Seif  %ustrahler%  eq  %maxstrahler%  Sethen 
ScS  maxstrahler  =  %maxstrahler%  +  1 
Self  %ustrahler%  gt  %maxstrahler%  Siithen 
Ses  maxstrahler  =  %ustrahler% 

Ses  ui  =  [calc  %ui%  +  1] 

Si:if  %ui%  le  %upnum%  Siithen;  &goto  innerloop  /******  Siiend 

/*  Silty  recoding  arc  %index%  %id% 
sel  %cover%-ID  =  %id% 
calc  shreve  =  %thisshreve% 
calc  strahler  =  %maxstrahler% 

Sclabel  continue  /*  finished  this  arc 
/*  Sety  that  was  arc  %index%  %id% 

&end  /*  finished  iteration  through  all  uncoded  arcs 
Si:goto  bigloop  /*  look  at  all  uncoded  arcs  again 
Siilabel  done 

save 
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Data  Listing 


Kr:Nr;if:Wim 
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stream  order  /  Basin  Relationship  study 


POLY_ 

BASIN_ 

NUM_ 

TOTAL 

TOTAL 

AVERAGE 

Observation  id  order  basin^rea 

PERIM 

STREAM 

LENGTH 

RELIEF 

SLOPE 

weM.-»P[;a 

Ytjriiywmn 
nwmmzti 

mmsma 

mtinmaa 
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SAS  Program 


OPTIONS  NOCENTER  PS=51  LS=80  NODATE  NONUMBER; 

FILENAME  streamdat  'C:\STREAMSTUDY\STREAM.DAT'; 

DATA  STREAMDATA; 

TITLE  'Stream  Order  /  Basin  Relationship  Study'; 

INFILE  streamdat  MISSOVER; 

INPUT  OBSIDPOLY_ID  ORDER  BASIN_AREA  BASIN__PERIM  NUM_STREAM  TOT_LENGTH 
TOT_RELIEF  AVG_SLOPE; 
logorder  =  log (order) ; 

RUN; 


***  DATA  LISTING  ****; 


RUN; 

PROC  PRINT  DATA=STREAMDATA; 
RUN; 


***  MULTIPLE  REGRESSION  ****; 

TITLE2  "MULTIPLE  REGRESSION  "; 

PROC  REG  DATA=STREAMDATA  LINEPRINTER; 

MODEL  ORDER  =  TOT_RELIEF  NUM_STREAM  TOT_LENGTH  AVG_SLOPE  /  influence  vif 
partial  collin  selection  =  backward; 

PLOT  RESIDUAL. *PREDICTED. ; 

OUTPUT  OUT=RESDATA  P=YHAT  R=E; 

RUN; 

***  RESIDUAL  ANALYSIS  ****; 

proc  univariate  data=RESDATA  normal  plot; 

var  e; 

run; 

QUIT; 
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REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  inforrnation  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  corrpieting  and  reviewing  this  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  Information,  including  suggestions  for 
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